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In this paper, we revisit the connection between the Riemann-Roch theorem and the zero energy 
solutions of the two-dimensional Dirac equation in the presence of a delta-function like magnetic 
field. Our main result is the resolution of a paradox - the fact that the Riemann-Roch theorem 
correctly predicts the number of zero energy solutions of the Dirac equation despite counting what 
seems to be the wrong type of functions. 
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INTRODUCTION 



The cross fertilization between mathematics and 
physics has proven itself fruitful over the history of sci- 
ence. Often important results in one field generate new 
ideas in others. For example, the seminal Atiyah-Singer 



index theorem in mathematics has spawned many novel 
directions of research in physics due to its connection 
with the Dirac equation with a background magnetic 
field. Examples include: string theory, particle theory, 
condensed matter theory, etc. While solving the Dirac 
equation with magnetic fields over arbitrary Riemann 
surfaces is an arduous task, the Riemann-Roch theo- 
rem, which counts functions of a particular type (namely, 
meromorphic functions), correctly predicts the analytic 
index of the Dirac Hamiltonian. When the magnetic field 
is sufficiently strong, this index agrees with the multiplic- 
ity of the zero energy solution. Therefore it is natural to 
suspect that the Riemann-Roch theorem can yield infor- 
mation for the E = eigenfunctions of the solutions to 
the Dirac equation. 

The celebrated Riemann-Roch theorem^ deals with 
functions that are analytic everywhere except a finite 
set of poles on a closed and compact Riemann surface 
X (a Riemann surface is a smooth, orientable surface). 
These functions of interest are known as meromorphic 
functions. The characteristics of meromorphic functions 
are summarized by a "divisor" D = ^2 peX n p Pi wnere 
{p} is a finite (and therefore discrete) set of points on 
X and n p are integers^. The vector space L(D) con- 
sisting of meromorphic functions that have poles of or- 
der < \n p \ at the points with n p < 0, and zeros of 
order > n p at for the points with n p > is the en- 
tity of interest. The Riemann-Roch theorem states that 
dim(L(D)) - dim(L(K - D)) = deg(D) + 1 - g with 
K any "canonical divisor" , g is the genus of X and 
deg(Z?) = J2 P ex n p> tne degree of the divisor. The con- 
cept of the canonical divisor is irrelevant for the purpose 
of this paper; this is because for divisors with sufficiently 
large degree, dim(L(K — D)) = 0, hence K does not en- 
ter. This is the situation we will be focusing on; however 
for further reading about canonical divisors, we refer the 
interested reader to Miranda or Hartshorne^. 

In physics it is known that the Riemann-Roch the- 
orem is connected to the analytical inde»2, of the two- 
dimensional Dirac equation under a constant magnetic 
field^. When the total number of flux quanta^ is suf- 
ficiently big, the latter reduces to the degeneracy of the 
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E = energy level. In this paper we show that this con- 
nection is quite subtle. First of all, meromorphic func- 
tions can diverge while the solutions 8 of the Dirac equa- 
tion in uniform magnetic field do not. This discrepancy 
can easily be removed by deforming the magnetic field so 
that 

B{f f )=J2 2m i s ^-n) (1) 

i 

while maintaining the total magnetic flux. After such 
deformation it is possible for the E = solution of the 
Dirac equation to diverge at a subset of {ri} 9 . The fact 
that such deformation preserves the analytic index of the 
Dirac operator is a consequence of the Atiyah-Singer in- 
dex theorem^. However a discrepancy remains even af- 
ter this deformation, namely, the solutions of the Dirac 
equation depends on both z and its conjugate z; these 
functions are not meromorphic. In this paper we ex- 
amine and resolve this discrepancy when X is the Rie- 
mann sphere. This case is particularly simple because for 
any divisor D with deg(D) > on the Riemann sphere, 
dim(L(A" — D)) = 0. Moreover, this is the a relevant 
example to realistic physical systems. Nevertheless, we 
believe that the idea exposed here can be generalized to 
other Riemann surfaces with higher genera. 



II. THE RIEMANN-ROCH PROBLEM ON THE 
SPHERE 

For the sake of completeness, in this section the mathe- 
matics of the Riemann-Roch problem will be introduced. 



A. Divisors 

In order to assist in the study of meromorphic func- 
tions, mathematicians introduce the concept of divisor. 

Definition II. 1. A divisor (Fig. Qp is a formal finite 
linear combination of points of X with integral coeffi- 
cients. The collection of divisors on a surface X, denoted 
Div(X) form an abelian group with the additive group 
operation defined as though distinct points were linearly 
independent vectors. Mathematicians define a partial or- 
dering on Div(X) by defining D > if all its coefficients 
are nonnegative. 

Divisors are of importance in the study of Riemann 
surfaces because the set of zeros and poles of a mero- 
morphic function on a surface X can be associated to a 
divisor. 

Definition II. 2. Let f be a nonzero meromorphic func- 
tion on a compact Riemann surface X. Then the the 
divisor of f is defined to be div(f) = J2 P ex n pP where 




FIG. 1: (color on-line) A graphical representation of a divisor. 
The white and black points denotes poles and zeros resp. 



n p is the order of the zero or pole of f at p. By conven- 
tion, n p is positive when f has a zero at p and negative 
when f has a pole at p. Such divisors are called principal 
divisors. 

It is important to note that such a divisor is always well 
defined because a meromorphic function on a compact 
Riemann surface can only have finitely many zeros and 
poles^. 

Definition II. 3. Let D be a divisor over a Riemann 
surface X, we define a complex vector space L(D) ~ 
{/ is meromorphic on X \div(f) + D > 0} 

The dimension of this space is the quantity of interest. 
While in general, it is very difficult to find the functions 
in i(-D), for the Riemann sphere, we can do so easily. 



B. L(D) and its Dimension on the Riemann Sphere 

First, we preform a stereographic projection to map 
the Riemann sphere onto Coo := C U {oo}, the one-point 
compactification of the complex planed. Without loss 
of generality, we may assume that in any divisor D, the 
coefficient associated to the point oo is 0. This can be 
done by appropriately choosing the point from which to 
stereographically project. 

Theorem II. 1. Let D = J2iLi n iPi £ Div(Coa) and 
deg(D) > 0. Letf D {z) = Y\ n l=l {z- Pi )-^ . Then L(D) = 
{g{ z )fD{z)\g{z) is a polynomial with deg(g) < deg(-D)} 

For those who arc not interested in the proof, read 
corollary III . 2 1 before proceeding to the following section. 

Proof. We will reproduce the proof in Miranda 1 page 149. 
Let g(z) be a polynomial with degree d. Then div(g) > 
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Corollary II.2. Let D £ Diw(C 00 ) with deg(D) > 0. 
Then dim(L(D)) = 1 + deg(-D). 



Proof. The set {f D (z)...,z n f D {z),...z de ^f D {z)} gen- 
erates L(D) as a complex vector space, therefore 
dim(L(L>)) = deg(L>) + 1 □ 

III. THE ZERO ENERGY SOLUTIONS OF THE 
DIRAC EQUATION IN A BACKGROUND 
MAGNETIC FIELD 



FIG. 2: (color on-line) The stereographic projection. The 
south pole is the base point. A point P on the sphere is 
projected to the plane by finding the intersection of the line 
connecting the north pole and P and the plane. 



-d oo, Hence we have 

div(f D ) = ^2 



-riiPi 



E 



Hi )O0. 



(2) 



We will now see the connection between the Riemann- 
Roch theorem and the Dirac equation. In order to do so, 
we first introduce the Dirac equation on two dimensional 
surfaces 3 . It is a well known result from topology that 
any two dimensional manifold can be constructed from 
a planar polygon with edges identified appropriately 11 . 
Therefore, writing the Dirac equation over these surfaces 
amounts to imposing additional periodic boundary condi- 
tions on the corresponding Dirac equation for flat space. 

We will use natural units by setting c = h = 1 as 
well as fixing the charge q of the particle so Dirac flux 
quantum 2irh/qc = 2ir. 



A. The Vector Potential and the Magnetic Field 



Therefore 

div(g(z)f D (z)) + D = div(g) + div(f D ) + D 

> TH - d)oo = (deg(D) - d)oo (3) 

i 

which is greater than or equal to if d < deg(D). 
Thus the space {g(z) f o{z)\g{z) is a polynomial with 
deg(g) < deg(Z?)} is a subspace of L{D). Now, for the 
other inclusion, take any nonzero h G L(D). Then 

div(h/fn) = div(h) — div{fu) 

>-D- div(f D ) = (-J2 n >)°° = - deg(£)oo (4) 



Let B(x,y) be a magnetic field. We define the vector 
potential A x (x,y) and A y (x,y) so that 



d x A y (x,y) - d y A x (x,y) = B(x,y) 



(5) 



Eq.© defines the vector potential up to a gauge trans- 
formation. The set of vector potentials which differ from 
each other by a gauge transformation form an equiva- 
lence class and may be considered the same. Choosing 
a representative of such an equivalence class is known as 
gauge fixing. In the following, we will fix the gauge by 
imposing 



which means that g can only have a pole of order at most 
deg(-D) at oo. Thus g is a polynomial of order at most 
deg(D). □ 

It is important to note that in the preceding equations, 
oo denotes the antipode of the point projecting to 0. This 
point arises because in general as \z\ — > oo, the function 
Jd{z) can have a pole or zero. 



d x A y (x,y) + d y A x (x,y) = 



(G) 



It will turn out that this choice of gauge will best exhibit 
interaction between the Dirac equation and complex ge- 
ometry. 

With the vector potential defined above, the time in- 
dependent Dirac equation in the presence of a magnetic 
field becomes^£ 
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-id x A x 4- d v — iA v 



o 



u(x,y) 
v{x,y) 



E 



u(x,y) 
v(x,y) 



(7) 




FIG. 3: (color online) As dictated by the Stokes theorem, the 
loop integral of the vector potential is equal to the magnetic 
flux. 



This integral equation translates into 



F(z)dz = 2iririi 



(12) 



Ci 



Now, by the Cauchy residue theoremi^, we can repre- 
sent F(z) as 



F(z) 



N 



(13) 



where Zi represents the position of the i-th magnetic flux 
point. The reader might have noticed that a holomorphic 
function G(z) could be added to F(z) while still preserv- 
ing the integral. Ignoring the holomorphic part of G{z) 
amounts to a further fixing of gauge. 



B. The Magnetic Flux and Divisors 



D. The Dirac Equation in Complex Coordinates 



The type of magnetic field given by Eq.([T]) admits a 
natural representation by a divisor, namely 



N 



D = ^2mpi<=> B(r) = 2TrJ2 n i S (r-r Pz ) (8) 



i=l 



C. The Vector Potential as a Complex Function 



Note that away from the r^s, B(x,y) = 0. Therefore, 
we have 



d x Ay dyA x — 
d x A x + dyAy = 



Let us define 



F(z, z) = A y (x, y) + iA x (x, y) 



(9) 



(10) 



Then Eq.© becomes d s F(z,z) = 0, the Cauchy- 
Riemann condition on F(z,z). As usual, d z — (d x — 
idy)/2, and d s = (d x + id y )/2. Moreover, by the Stokes 
theorem, the equation d x A y — d y A x = 27r£\ rii5(r — fi) 
rewrites into the integral form 



A(x, y) ■ df = 2-nrii 



(11) 



Ci 



where Ci is a counter-clockwise loop enclosing only the 
i-th magnetic flux point (Fig. [3]) 



We first make the following standard definitions: 
1 



A, 



{A x — iAy 



A z — 



(14) 



Using z and z as coordinates, the Dirac equation in the 
presence of a magnetic field reads 
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dz - iA* 



d z 



iA 2 







u(z, z) 
v(z,z) 



E 



u(z, z) 
v(z,z) 



(15) 



Taking E = 0, Eq.^ becomes: 

(d s - iA z )u(z,z) = 
(8 z -iA z )v(z,z)=0 (16) 

Substituting in Eq. (|10[) . these two equations become 



(d z + -F(z))u(z,z) = 
(d,-±F(z))v(z,z) = 



(17) 



E. Solving the Dirac Equation 

The second line of Eq. © ensures that it is possible to 
pick a &(x,y) so that 



A : , 



-dy$ Ay=d x $, 



(18) 
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or equivalently, A z = —id z Q>/2 and A z = id z &/2. Since 
A z = —iF/2 and A z — iF/2, this means 



1 N 



1 



2 ^ z - zi 2 

I— 1 



2 ^ 2 - zi 2 

i—1 



(19) 



In writing down the above equations we have used 
Eq.CJJJ. 



(20) 



A simple calculation shows that 

N 



i=l 



satisfies the above equations. Substituting $ into the 
Dirac equation allows us to write the two differential 
equations as 



(d M + d M $(z,z))u(z,z) = 
(d z -d z Mz,z))v(z,z) =0 



(21) 



Solving this equation by separation of variables gives the 
solutions 



u(z,z) = f(z)e-*W 
v(z,z)=g(z)e*M 



(22) 



Notice that the functions f{z) and g(z) do not contribute 
to differentiation with respect to z and z respectively. 
Substituting in <J>(z, z) gives 



x(z,z)= / (,)n^- F 

N 



(23) 



Because our original divisor does not specify any poles 
or zeroes at oo, v(z, z) is not a valid solution because 
it diverges as \z\ — > oo. Now in order for u to remain 
bounded at infinity, f(z) must be a polynomial satisfy- 
ing deg(/(z)) < The reason that f(z) cannot be 
a rational function of the form, say l/z m , is because it 
makes the divisor smaller under the partial ordering of 
Div(X) (see "Divisors"). Thus we conclude that there 
are 1 + rii linearly independent solutions obtained 
by setting f(z) equal to 1, z, z 2 ...z 1+ ^i ni . Notice that 
1 + Y^i n i = deg(.D) + 1, a number which we calculated 
to be the complex dimension oiL{D) from a purely math- 
ematical viewpoint in the previous sections. 

Despite the above success, the function u{z, z) in 
Eq. (|23[) is not meromorphic hence is not a member of 
L(D). Another equivalent statement is that Eq. pTj) is 
not the Cauchy-Riemann equation d z u = 0. We will 
remedy this situation with the use of a phase transfor- 
mation. 



IV. THE PHASE TRANSFORMATION 



Our objective now is to transform functions u(z, z) = 
f( z ) YliLi | z _^.|n i into meromorphic functions. This can 
be done by multiplying u by the unit modulus complex 
valued function x(z, z) = Ililid^f 1 )"^ 2 - B y taking the 
product x u i we S e t the meromorphic function 

N 1 

X (z,z)u(z,z) = f(z)l[- — . (24) 

Where f(z) = 1, z, ni . Notice that this gives ex- 
actly the space L(D) of Theorem 2.1. 



A. Monodromy 



We now show \ is single valued. It is in fact a nontrivial 
result that \ does not have any branches. Because \ nas 
unit modulus, we can write it as e l6 ^ z ' z K After doing so, 
it suffices to study the behavior of this function along a 
loop around any one of the Zk ■ 

Theorem IV. 1. As x(z,z) wraps once around any one 
of the Zk, 9{z 1 z) changes by —2imk (an example is shown 
in Fig^j. 



Proof. First it is clear that if the loop encloses only 
Zk, the only factor in \ that will exhibit winding is 
W 2 . Now, we have 



z - z k = \z - z k \e 
z- Zk=\z- z k \e 

Therefore, we have 



iarg(z-z fc ) 
-iarg(z-z fc ) 



Z - Z k 

z- z k 



_ e ~in k arg(2-z fc ) 



(25) 



(26) 



Now clearly as z wraps around Zk once counterclockwise, 
the exponent of the right hand side changes by —2i:nk 

□ 



Now we can demonstrate that this phase transforma- 
tion is single valued. 

Corollary IV. 2. ~x(z, z) has no branches. 



Proof. After a full loop A6 = -2-Kn k , so e lAe = 1 thus 
after a loop, the value of x(z, z) remains the same. □ 
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FIG. 4: (color on-line) A plot of the real and imaginary parts 
of x{ z i z) as a vector field. The associated with the red/blue 
points are ±1. 



B. Recovery of Cauchy-Riemann Equations from 
the Dirac Equation 



We will now show that the transformation of u(z, z) 
into a meromorphic function is not a coincidence. It 
results from the fact that the Dirac equation trans- 
forms into the Cauchy-Riemann equations after the phase 
transformation. 

Theorem IV. 3. The equation (dg + dz$)u = trans- 
forms into dz{x u ) — 0- 

Proof. Let ((z) = \{ z , z)u(z , z) . Then u — x~ 1 C — xC- 
We have 



d- z u = dz{xQ = (dzx + xdsC 



(27) 



By an elementary computation x^zX — — dz*&- As a re- 
sult, the Dirac equation becomes 



(dzx = CMx = -dM(x) 



(29) 



(28) 



This cancels the extraneous dz&((x) m t- ne previous 
equation giving us the Cauchy-Riemann equation. □ 



CONCLUSION 



On the Riemann sphere, we have shown that upon 
a phase transformation the zero energy solutions of the 
Dirac equation in the presence of a magnetic field given 
by Eq.|T]) becomes the meromorphic functions of the 
Riemann-Roch theorem. We have also performed the 
analogous calculation for the complex torus in which 
there is a periodic boundary condition requiring us to 
represent the vector potential in terms of the elliptic 
theta function 1 . In this case, the phase transformation 
is replaced by the ratio of elliptic theta functions. 

In general, for Riemann surfaces of higher genera, it 
is not clear that such phase transformations exist. This 
is because the phase transformation must be consistent 
with the appropriate periodic boundary conditions dis- 
cussed earlier. A related fact is that it is difficult to 
prove the existence of nonconstant meromorphic func- 
tions on arbitrary compact Riemann surfaces; this was 
first demonstrated by Riemann^. We conjecture that 
the existence of a phase transformation is equivalent 
to the existence of nonconstant meromorphic functions. 
If this is true, there might be an alternative proof for 
the fact that compact Riemann surfaces are "projective 
algebraic"—. 
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